Control of defect-mediated tunneling barrier heights in ultrathin MgO films 
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The impact of oxygen vacancies on local tunneling properties across rf-sputtered MgO thin films was inves- 
tigated by optical absorption spectroscopy and conducting atomic force microscopy. Adding O2 to the Ar 
plasma during MgO growth alters the oxygen defect populations, leading to improved local tunneling char- 
acteristics such as a lower density of current hotspots and a lower tunnel current amplitude. We discuss a 
defect-based potential landscape across ultrathin MgO barriers. 



PACS numbers: 81.15.-z, 61.72.jd, 85.30.Mn 

MgO-based tunnel barriers have been intensively stud- 
ied since the theoretical prediction^] of >1000% tun- 
neling magnetoresistance (TMR) in Fe/MgO/Fe and ini- 
tial experimental confirmation. 3 4 However, TMR ratios 
at room temperature in MgO-based MTJs have reached 
only 200% with Fe electrode^ and 600% with Fe- 
CoB alloy electrodes.^ The most important cause of the 
discrepancy is believed to be structural defects in the 
MgO barrier. Yet despite previous experimental"^ and 
theoretical^^ studies of point defects in MgO, the prop- 
erties of defects in MgO films are still in debate. As 
technology goes nanoscale and increasingly integrates al- 
ternative dielectrics to silicon, it becomes important to 
understand the impact of defects on the properties of 
these ultrathin dielectrics such as MgO toward coherent 
transport 17 and memristive effects EMEU 

To alter the density of only intrinsic MgO defects in 
our films, 21 we fixed the working pressure during rf- 
sputtering but modified the proportion of O2 and Ar from 
0% to 10% O2. To prevent the oxidation of the Fe sur- 
face, the first half monolayer was sputtered with pure 
Ar gas for all samples. 20 To qualitatively determine the 
change in defect density, we measured the optical trans- 
mission of MgO (001) / /MgO (50 nm) samples using a UV- 
Vis-NIR Perkin-Elmer Lambda 950 spectrophotometer 
with an integrating sphere of 150 mm. Local tunneling 
current maps and cur rent- volt age (I-V) curves were ob- 
tained on MgO(001)//MgO(10 nm) /Fe(20 nm) /Mg O (1.2 
nm) samples grown in optimized growth condition d 21 * 22 * 
using conductive atomic force microscopy (C-AFM; N- 
Tracer, Nanofocus Inc.). The C-AFM tip was grounded 
and the lower Fe electrode was biased. The spatial res- 
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olution of C-AFM is the same as that of a conventional 
AFM and the minimum detectable current of C-AFM is 
about 0.5 pA. 

The spectral dependencies of the optical absorption co- 
efficient a(uj) of MgO(OOl)/ /MgO(50 nm) samples shown 
in Fig. 1 were extracted from experimental transmission 
spectra using ai ayer (uj) = - [In T tota i(uj) - In T sub (uj)]/d, 
where T(u) is the optical transmission and d is the MgO 
layer thickness atop the MgO substrate with its own opti- 
cal transmission T sv b(uS). Note the presence in the sam- 
ple with 0% O2 of three peaks (see arrows in Fig. 1): 
the peak at ~4.2 eV is assigned to M centers (the com- 
bination of two oxygen vacancies); 11 the peak at ~5.0 
eV is assigned to F and F + cen ters (neutral and sing ly 
charged oxygen vacancies) GEH3I finally, a weak peakP^at 
~5.6 eV i s ass igned to F 2+ center (doubly charged oxygen 
vacancy) P^U As the O2 proportion increases, the spectral 
weight around each of the three peaks decreases system- 
atically, underscoring a control during sputter-growth of 
oxygen defect densities in our MgO films that is more 
direct and effective than other methods fHHHIHl 
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FIG. 1. (Color online) Optical absorption spectra of 
MgO(001)//MgO(50 nm) samples. 



2 



V 




FIG. 2. (Color online) (a) Topographic map of MgO barrier 
in 0%-O2 MgO layer and local tunneling maps with 1 V bias 
of (b) 0%- and (c) 10%-O 2 MgO layers. Topographic map 
in 10%-O2 MgO layer is nearly the same with (a). The rms 
roughness of (a) is 0.14 nm. 



We now examine the impact of oxygen vacancies on 
tunneling transport across a 1.2 nm-thick MgO layer. 
Figures 2(a) and 2(b) respectively show the topographic 
and local tunnel current maps across a 0%-O2 MgO layer 
with 1 V bias. The MgO layer is very flat owing to a rms 
roughness that is identical to th at of the MgO substrate. 
Similar to previous results pfflMUl we obtain an inhomo- 
geneous local current map including a number of current 
hotspots. Since there is no correlation between topogra- 
phy and current maps, this tunnel current inhomogeneity 
originates from fluctuations in barrier thickness and/or 
height. 28 However, while the topographical map of the 
10%-O 2 MgO layer resembles that of the 0%-O 2 MgO 
layer (data not shown), the corresponding local tunnel 
current map [see Fig. 2(c)] reveals a lower density of 
hotspots and, moreover, a lower tunnel current ampli- 
tude. This suggests that the tunnel current inhomogene- 
ity across MgO reflects solely a barrier height fluctuation 



that originates from oxygen vacancies. 

In support of this suggestion, we show in Fig. 3(a) 
statistical distributions of local tunnel current through 
0%- and 10%-O 2 MgO ultrathin films. If there is a 
distribution of tunneling parameters, such as barrier 
height /thickness, the local tunnel current is correctly de- 
scribed by a log-normal distribution. The tunnel current 
i = i§e~ l l x depends exponentially on the barrier thick- 
ness I and the attenuation length A, where A is a func- 
tion of barrier height. Assuming a Gaussian distribution 
of l/X with a mean \i and a standard deviation a. th e 
tunnel current has a log-normal distribution P(z)p3SH 
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By fitting our data for 0%- and 10%-O 2 MgO layers [see 
lines in Fig. 3(a)], we respectively find typical currents 
i typ of 0.38 pA and 0.047 pA, averaged currents (i) of 6.1 
pA and 0.30 pA, and the a values of l/X of 1.4 and 1.1. 
Since both MgO layers have virtually the same thickness 
fluctuation [see Fig. 3(b)], a lower a value implies a spa- 
tially more uniform barrier height, and thus the lower i typ 
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FIG. 3. (Color online) (a) Distributions of local tunnel cur- 
rent from Figs. 2(b) and (c). Lines are the best fits with 
a log-normal distribution, (b) Distributions of topographic 
deviation from the nominal thickness of MgO barriers in 0%- 
and 10%-O2 MgO layers. These distributions fit well with a 
Gaussian function, where the standard deviation is 0.26 and 
0.25 for 0%- and 10%-O 2 MgO layers, respectively. 
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FIG. 4. (Color online) Current- voltage curves of 0%- and 
10%-O2 MgO layers, measured at hot spot and background 
using conductive tip. 



and (i) imply a larger barrier height as the O2 proportion 
is increased. 

To understand this increased barrier height from a de- 
fect perspective, we present in Fig. 4 typical local I-V 
curves taken away from and on current hotspots. 29 Since 
the first MgO half monolayer was systematically grown in 
a pure Ar gas to avoid Fe oxidation, the I-Vs for the 10%- 
O2 MgO may be asymmetric, with data at negative bias 
(thus probing the bottom interface) that resembles that 
for a 0%-O2 MgO sample. We therefore focus only on 
data acquired at positive bias, corresponding to the top 
interface, and qualitatively pinpoint the barrier heights. 
Going from 0% to 10% O2 raises the barrier height from 
0.9 eV to 2.3 ~ 4.1 eV in the background, 30 and from 
0.4 eV to 1.1 eV on hotspots. Interestingly, the I-Vs for 
0%-O2 /'off- hotspot' and 10%-O2/'on- hotspot' are quite 
similar. This suggests that going from 0% to 10% O2 
eliminates hotspot-inducing defects responsible for the 
0.4 eV barrier height in the 0%-O2 sample, while the de- 
fects responsible for the 0.9 eV background barrier in the 
0%-O2 sample now account in the 10%-O2 sample only 
for hotspots relative to an improved background. In fact, 
for a given positive bias, the tunneling current is lower for 
10%-O2/'on-hotspot' than for 0%-O2/'off-hotspot', im- 
plying that the density of the defect type(s) responsible 
for this I-V characteristic has decreased. Naively, if we 
assume Fermi level pinning in the middle of the bulk MgO 
band gapl 31 | S2 1 of 7.8 eV and the presence at 5 eV above 
the top of the valence band of F/F + defect states^ then 
F/F + defects define an effective barrier height of ~1.0 eV 
that might correspond to that observed in the I-V data 
for the 0%-O2 background and 10%-O2 hotspot data, 
consistently with the near extinction of F/F + absorption 
at 5 eV as the oxygen proportion is increased (see Fig. 1). 
A more quantitative description of the defect-mediated 
potential landscape across ultrathin MgO would require 
a) mitigating the impact of water exposure between the 
growth and characterization phases of the study; b) de- 



termining the evolution of the band gap across the atomic 
planes of the ultrathin film so as to account for metal- 
induced gap states^ and size effects; 33 and c) considering 
the likely surface/interface nature of defects in the ultra- 
thin film, whose transition energies and positions within 
the band gap are different from those of the bulk and 
remain largely unstudied. 11 

In conclusion, we have shown how to control the na- 
ture and density of oxygen defects during the growth of 
rf-sputtered MgO thin and ultrathin films via the combi- 
nation of optical absorption and local tunneling current 
maps obtained by conductive atomic force microscopy. 
Eliminating certain defect populations leads to improved 
barrier heights with a reduced spatial fluctuation, and 
to a change in the type of defect-induced hotspot that 
drives tunneling transport in a micronic device. 

We thank Y. Yin for assistance with optical exper- 
iments, and D. Halley and S. J. Park for stimulating 
discussions. This research was supported by the EC 
FP6 (NMP3-CT-2006-033370), the French ANR (ANR- 
06-NANO-033-01, ANR-09-JCJC-0137), and by the Ko- 
rean MEST (2010K000339 from CNMT under 21st Cen- 
tury Frontier R&D Programs; NRF 2010-0020416). 

1 W. H. Butler, X.-G. Zhang, T. C. Schulthess, and J. M. Maclaren 

Phys. Rev. B 63, 054416 (2001). 
2 J. Mathon and A. Umerski, Phys. Rev. B 63, 220403R (2001). 
3 M. Bowen et al, Appl. Phys. Lett. 79, 1655 (2001). 
4 J. Faure- Vincent et al, Appl. Phys. Lett. 82, 4507 (2003). 
5 S. Yuasa, T. Nagahama, A. Fukushima, Y. Suzuki, and K. Ando, 

Nature Mater. 3, 868 (2004). 
6 S. Yuasa, J. Phys. Soc. Jpn. 77, 031001 (2008). 
7 S. Ikeda et al, Appl. Phys. Lett. 93, 082508 (2008). 
8 Q. S. Wang and N. A. Holzwarth, Phys. Rev. B 41, 3211 (1990). 
9 A. Gibson, R. Haydock, and J. P. LaFemina, Phys. Rev. B 50, 

2582 (1994). 

10 F. Illas and G. Pacchioni, J. Chem. Phys. 108, 7835 (1998). 
11 D. Domingues-Ariza, C. Sousa, F. Illas, D. Ricci, and G. Pac- 
chioni, Phys. Rev. B 68, 054101 (2003). 
12 J. P. Molnar and C. D. Hartman, Phys. Rev. 79, 1015 (1950). 
13 H. R. Day, Phys. Rev. 91, 822 (1953). 

14 Y. Chen, T. Williams, and W. A. Sibley, Phys. Rev. 182, 960 
(1969). 

15 P. G. Mather, J. C. Read, and R. A. Buhrman, Phys. Rev. B 73, 
205412 (2006). 

16 G. H. Rosenblatt, M. W. Rowe, G. P. Williams, Jr., R. T. 
Williams, and Y. Chen, Phys. Rev. B 39, 10309 (1989). 

17 F. Greullet et a/., Phys. Rev. Lett. 99, 187202, (2007). 

18 M. Bowen et al, Appl. Phys. Lett. 89, 103517 (2006). 

19 D. Halley et al, Appl. Phys. Lett. 92, 212115 (2008). 

20 Any Fe oxidation would lead to states near the middle of the 
MgO band gap, and thus to a lower barrier height in contrast 
to our results. See P. A. Cox, Transition Metal Oxides, Oxford 
Science Publications, Oxford, Chap. 2 (1992). 

21 D. J. Kim et al, J. Phys. D: Appl. Phys. 43, 215003 (2010). 

22 Careful thickness calibrations using a combination of ex-situ sty- 
lus profilometry and an in-situ quartz microbalance enable an 
effective relative statistical error of 0.1 nm on our 1.2 nm-thick 
films. 

23 Low resolution reflects here a loss of transmission across MgO. 
24 G. X. Miao, Y. J. Park, J. S. Moodera, M. Seibt, G. Eilers, and 

M. Munzenberg, Phys. Rev. Lett. 100, 246803 (2008). 
25 V. Da Costa, Y. Henry, F. Bardou, M. Romeo, and K. Ounadjela, 

Eur. Phys. J. B 13, 297 (2000). 
26 K. M. Bhutta, J. Schmlhorst, and G. Reiss, J. Magn. Magn. 

Mater. 321, 3384 (2009). 



4 



Y. Ando, H. Kameda, H. Kubota, and T.Miyazaki, J. Appl. Phys 
87, 5206 (2000). 

F. Bardou, Europhys. Lett. 39, 239 (1997). 

'Nonlinear data excludes direct conduction between the tip and 
the lower Fe layer even on the hotspot. 

'Depending on the local background spot considered, current in- 
creases may appear at 2.3, 3.1, 3.6, and 4.1 eV. 
M. Bowen et al, Phys. Rev. B 73, 140408R (2006). 



This is consistent with barrier heights of 3.6 eV and 3.9 eV found 
for MgO/Fe and MgO/Au. See W. Wulfhekel et al, Appl. Phys. 
Lett. 78, 509 (2001) & S. Guezo et al, Appl. Phys. Lett. 93, 
172116 (2008). 

A. Ohtomo, D. A. Miiller, J. L. Grazul, and H. Y. Hwang, Nature 
419, 378 (2002). 



